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ABSTRACT 

A substantial fraction of barred galaxies host additional nuclear bars which tumble with 
pattern speeds exceeding those of the large-scale (primary) stellar bars. We have investigated 
the mechanism of formation and dynamical decoupling in such nested bars which include gaseous 
(secondary) nuclear bars within the full size galactic disks, hosting a double inner Lindblad 
resonance. Becoming increasingly massive and self-gravitating, the nuclear bars lose internal 
(circulation) angular momentum to the primary bars and increase their strength. Developing 
chaos within these bars triggers a rapid gas collapse — bar contraction. During this time period, 
the secondary bar pattern speed fl s ~ a -1 , where a stands for the bar size. As a result, fl s 
increases dramatically until a new equilibrium is reached (if at all) , while the gas specific angular 
momentum decreases — demonstrating the dynamical decoupling of nested bars. Viscosity, and 
therefore the gas presence, appears to be a necessary condition for the prograde decoupling of 
nested bars. This process maintains an inflow rate of ~ 1 M yr _1 over ~ 10 s yrs across the 
central 200 pc and has important implications for fueling the nuclear starbursts and AGN. 

Subject headings: hydrodynamics — galaxies: active — galaxies: evolution — galaxies: kinematics and 
dynamics — galaxies: spiral — galaxies: starburst 



1. Introduction 

A large fraction, ~ 1/3, of barred galaxies host 
a secondary (nuclear) bar in the central regions, 
in addition to the primary, large-scale stellar bar 
(Laine et al. 2002; Erwin & Sparke 2002). Iron- 
ically, such nested bar systems have been pre- 
dicted theoretically, including their most intrigu- 
ing property — the state of dynamical decoupling 
(Shlosman, Frank & Begelman 1989; Shlosman, 
Begelman & Frank 1990). Few examples of dou- 
ble barred galaxies have been observed already in 
the 1970s, but have been classified as having tri- 
axial (elliptical) bulges rather than nuclear bars 
(de Vaucouleurs 1974; Sandage & Brucato 1979; 



Kormendy 1982). Both stellar and mixed type 
bars (e.g., Buta & Crocker 1993; Shaw et al. 
1995; Friedli et al. 1996; Jungwiert, Combes & 
Axon 1997; Jogee, Kenney & Smith 1998; Knapen, 
Shlosman & Pclcticr 2000), as well as gaseous nu- 
clear bars (e.g., Ishizuki et al. 1990; Devereux, 
Kcnncy & Young 1992; Forbes et al. 1994; Mirabel 
et al. 1999; Kotilainen et al. 2000; Maiolino et al. 
2000) have been detected. Random mutual orien- 
tation of the observed nested bars confirms that 
they spend a substantial fraction of their lifetime 
in the decoupled state — tumbling with different 
pattern speeds. 

Numerical simulations have shown indeed that 
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pattern speeds of the mixed nuclear bars are sub- 
stantially larger than those of the primary stel- 
lar bars (e.g., Friedli & Martinet 1993: Combes 
1994; Heller & Shlosman 1994). However, purely 
stellar nuclear bars form only under special initial 
conditions (e.g., Friedli & Martinet 1993), while 
gaseous bars have been claimed not to decouple 
(Wada & Habe 1992; Combes 1994). A number 
of issues have never been resolved — what trig- 
gers the dynamical decoupling of nested bars? Do 
purely gaseous self-gravitating bars form, and can 
they decouple from the primary bars? In this Let- 
ter we use numerical simulations of self-gravitating 
gas in galactic disks to show that gaseous bars 
form and decouple as a result of stellar bar-driven 
inflow into the central kpc. Furthermore, we study 
the physical reasons for this dynamical runaway in 
nested bar systems (see also Shlosman 2003; and 
in preparation). 

2. Numerical Modeling 

Cold gas forms a thin layer in disk galaxies. Be- 
cause we are basically interested in the dynamical 
evolution of this layer, we limit the simulations to 
the disk midplane and use the updated 2-D version 
of our Smoothed Particle Hydrodynamics (SPH) 
code (e.g., Heller & Shlosman 1994). For test pur- 
poses we also make use of the grid code ZEUS-2D 
(Stone & Norman 1992). The initial conditions 
have been tailored in order to focus on the role of 
central resonances in disk galaxies by modifying 
the background analytical potential of Miyamoto 
& Nagai (1975), $(r) = -GM/Vr 2 + A 2 , using 
prescription from Englmaier & Shlosman (2000). 
Here M is the mass of a halo, bulge or a disk, and 
A is the radial lengthscalc, all given in Table 1, 
and r is the radius vector in the disk plane. Ad- 
justments to this potential allow us to specify the 
positions of axisymmetric resonances: the inner 
inner and outer inner Lindblad resonances (IILR, 
OILR) and the corotation. Positions of these res- 



Table 1: Model Parameters 



Component 


Mass (10 1U M ) 


A (kpc) 


halo 


5 


9 


bulge 


0.4 


0.2 


disk 


2 


3 


gas 


0.05 


3 



onances correspond to the solutions of equation 
fl — k/2 = fl p , where O, k and O p are angular ve- 
locity, epicyclic frequency and primary bar pattern 
speed. For the models presented here, those are 
fixed at r = 0.5 kpc, 2 kpc and 6 kpc respectively. 
The primary stellar bar is given by the Ferrers' 
(1877) potential with n = 1 and rotates with a 
prescribed pattern speed f2 p = 25.46 Gyr -1 . Its 
mass, the semi-major (a) and semi-minor (b) axes 
are 4 x 10 9 M©, 5 kpc and 1.25 kpc, respectively. 
The axisymmetric contribution from the primary 
bar mass is included in the initial conditions. The 
non-axisymmctric part is brought up gradually in 
order to avoid transients. The number of SPH 
particles has been varied between N = 10 4 and 

4 x 10 5 and the minimal gravitational softening 
has been fixed at 150 pc. A large number of mod- 
els have been run with increasing gas contents. 
Here we only discuss a representative model show- 
ing a nested bar descoupling, with the gas mass 

5 x 10 s M Q within the corotation, which corre- 
sponds to <~ 2.7% of the total mass within 3 kpc 

- the initial radial scalclcngth of the isothermal 
gas. The gas sound speed is 10 km s _1 . 

3. Results 

As the large-scale bar is brought up gradually, 
the gas within the corotation radius forms a pair 
of trailing offset shocks and falls inwards losing 
its angular momentum via gravitational torques 
induced by the primary bar (see Fig. 1 and An- 
imation Sequence 1). At around ~ 1 kpc, the 
shocks curve azimuthally, and the gas flow is 
shock-focused into more circular orbits attempt- 
ing to form an ovally-shaped nuclear ring. The 
gas inside the central kpc is losing its angular mo- 
mentum forming a secondary bar which leads the 
primary bar by a constant angle ~ 50° — 60° and a 
pattern speed O s = fi p , and has a flat density dis- 
tribution. A pair of leading shocks form between 
this nuclear bar and the ring. Almost immediately 
the ring fragments and gradually merges with the 
nuclear bar. At around t ~ 0.35 — 0.4 Gyr this 
secondary bar with a semimajor axis of 1 kpc and 
ellipticity e ~ 0.65 decouples in the retrograde di- 
rection, in the primary bar frame (still prograde 
in the inertial frame). One can follow up the sub- 
sequent evolution in one of the three frames of ref- 
erence: inertial, primary or secondary. We mostly 
employ the first two. 
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Fig. 1. — Gas evolution in the central 3 kpcx3 kpc part of a barred galaxy in the frame of the primary bar 
(Animation Seq. 1). Rotation is anti-clockwise and the primary bar is horizontal. The time in Myrs is shown 
in the right upper corners (see text for more details). 



This unusual behavior of the gaseous bar is a 
result of gravitational torques from the primary 
bar (Heller, Shlosman & Englmaier 2001; see also 
below). Its swing towards the position angle of 
the primary bar is accompanied by a loss of in- 
ternal (i.e., circulation) and precession (tumbling) 
angular momenta in the gas. The bar ellipticity 
increases to <~ 0.85 (Fig. 2) and f2 s drops to <~ 
12 Gyr -1 (Fig. 3) in the incrtial frame. Next, the 
gaseous bar reverses its tumbling to the prograde 
one (in the primary frame) and quickly shrinks to 
much smaller radii. This reversal coincides with 
a dramatic gas inflow (Fig. 4). The gas density 
distribution in the secondary bar becomes cen- 
trally peaked. The avalanche inflow clearly leads 
to a period of an unstable bar shape, <~ 0.1 Gyr 
long. The bar apparently fragments during the 
collapse and two fragments preserve their iden- 
tity and store some of the circulation angular mo- 
mentum during this time (e.g., t — 0.474 Gyr 
in Fig. 1). The final value of Q s ~ 85 Gyr -1 is 
achieved at t ~ 0.6 Gyr with some superposed bar 
length oscillations which are well correlated with 



the mutual orientation of both bars. 
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Fig. 2. — Evolution of ellipticity, e = 1 — b/a (dot- 
ted line), and semimajor axis, a (solid line) of the 
secondary bar. 
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Fig. 3. — Pattern speeds of primary (f2 p , dotted 
line) and secondary (fi s , solid line) bars in the 
inertial frame. 
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Fig. 4. — Evolution of gas inflow rates across the 
central 175 pc (solid line), 600 pc (dashed line) 
and 1 kpc (dotted line). The peak at 0.45 Gyr is 
unresolved. 

The observed pattern speed reversal (i.e., 
speedup in the inertial frame) represents the dy- 
namical runaway (decoupling) of the gaseous bar 
from the primary bar. The secondary pattern 



speed has first decreased from f2 s = f2 p to its min- 
imal value of ~ 12 Gyr -1 and then increased to the 
maximal <~ 85 Gyr -1 ~ 3.4il p where it has stabi- 
lized (Fig. 3). By t <~ 0.7 Gyr the bar has shrank 
to ~ 150 pc. This final bar size is the result of a 
limiting gravitational softening (Section 2). Fig. 4 
shows the gas mass inflow rates across a number 
of characteristic radii. The time delays between 
the peak rates confirm that the avalanche inflow 
propagates from larger to smaller spatial scales. 
The mass ratio of secondary-to-primary bar, /, 
toward the end of this evolution is ~ 0.1, and the 
gaseous bar contains about 16% of the total mass 
within its radius. 

Next, we analyze this gas evolution and address 
the focal questions — what exactly triggers the 
secondary bar decoupling process and what deter- 
mines the final secondary pattern speed. We only 
focus on the essential details here. Some of the 
observed features in this model do not seem to 
be of a principal importance. For example, the 
retrograde swing of the gaseous bar, although is 
interesting and maybe relevant in some systems, 
is not crucial for the subsequent prograde decou- 
pling as other models show (in preparation). On 
the other hand, we find that changes in the gas 
density and bar size lead directly to the increase 
in the central mass concentration in the bar, thus 
amplifying the self-gravitational effects in the gas, 
affecting dynamics of the system, and triggering 
the observed runaway. 

First, when the secondary bar forms, it ex- 
tends to the radius of the Max(f2 — k/2) curve, 
at <~ 1 kpc. This is also the radius where the nu- 
clear ring attempts to form. The radial density 
profile in the gaseous bar is flat and even centrally 
depressed, because the IILR not only prevents the 
gas inflow to within 500 pc of the center, but also 
applies torques which expcll some of the gas across 
the IILR. The gas response between the ILRs is 
<~ 50 — 60° out of phase with the primary bar 1 
defining the initial position angle of the secondary 
bar. The Jacobi energies, Ej (e.g., Binney & 
Tremaine 1987), of the gas are distributed around 
the IILR as a Gaussian with AEj /Ej <~ 5% - 10% 

The so-called periodic X2 orbits in the notation of Con- 
topoulos & Papayannopoulos (1980) arc oriented at 90° to 
the bar and serve as attractors for the gas motion. The 
viscous torques prevent the gas response from being com- 
pletely aligned with the X2 orbits 
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(e.g., Fig. 4 of Heller et al. 2001). 

Second, two important developments accom- 
pany this evolution: the gaseous bar mass in- 
creases due to the continuous inflow, and the spe- 
cific internal angular momentum (i.e., circulation) 
in the bar decreases due to the ongoing torquing 
from the primary bar (note that the secondary 
bar is positioned in the first and third quadrants) . 
Loss of internal circulation increases the bar cllip- 
ticity by diminishing its minor axis, thus forcing 
the gas inwards across the IILR (i.e., as measured 
by Ej distribution of gas particles). When <~ 50% 
of the gas mass crosses the IILR, it does not find 
any supporting periodic orbits — these are aligned 
with the primary bar inside the IILR. The primary 
bar then torques the gas, which, because it leads 
the bar by less than 90°, swings in the retrograde 
fashion towards the main bar. Fig. 3 clearly dis- 
plays this Sl s slowdown in the inertial frame. 

Third, increasing e (i.e., bar strength) to 0.8 
in the stellar bars generates a substantial chaos 
(e.g., Udry & Pfenniger 1988; Martinez- Valpuesta 
& Shlosman 2004) — the self-gravity is necessary 
for this. For a gaseous bar, this leads to a shock 
dissipation and to a dramatic gas inflow. El-Zant 
& Shlosman (2003) have studied the stability of 
nested bar systems using the Liapunov exponents 
and find that this system has a narrow "window" 
for a stable co-existence — otherwise a mutually 
excited chaos develops. If / <S 0.01, the secondary 
bars cannot exist in the decoupled phase as no 
orbits aligned with these bars have been found. 
On the other hand, when / <; 0.15, the primary 
bar becomes increasingly chaotic at the bar-bar 
interface. The value / ~ 0.1 obtained here fits 
comfortably within the allowed range. 

The main effect of intrinsic dissipation is the 
sudden increase in the central mass concentration 
and in the gas self-gravity. This is nicely exhibited 
by the radius of Max(f2 — n/2) which is dragged in- 
wards by the collapsing gas from ~ 1 kpc. In fact, 
the secondary bar length follows exactly this ra- 
dius, which represents the position of the ILRs to a 
large degree. This is not surprising, as these reso- 
nances are main generators of chaos in disk galax- 
ies. Because the gas cannot reside on intersecting 
orbits, gaseous bars are expected to be limited in 
size to within their ILRs, if the latter exist. 

Although the gas in the model contributes only 
<~ 2% of the total mass within the primary bar 



corotation at 6 kpc, it alters substantially the po- 
sition of the IILR at later times, as the gas fraction 
of the total mass within the secondary bar radius 
ends up at ~ 16%. This is verified through the ap- 
pearance of a characteristic orbital family in the 
nonlinear orbit diagram (in preparation), but even 
the linear (epicyclic) diagrams exhibit this clearly. 

Overall, the gaseous bar shrinks by a factor of 
<~ 7 from the decoupling moment and its pattern 
speed fi s increases by the same amount (Figs. 2, 
3). Their product, fl s a ~ const., varies by ±22% 
only. This corresponds to a decrease in the tum- 
bling angular momentum of the bar, ~ fl s a 2 , by 
a factor of <~ 7. Qualitatively this means that 
the gaseous bar is speeding up when contracting, 
and is losing some of its angular momentum to the 
primary bar torques. To verify that this effect is 
driven by the gas self-gravity, we have evolved the 
same model neglecting the gas gravity. The evolu- 
tion diverges at the retrograde swing — the non- 
sclf-gravitating bars do not generate additional 
chaos and dissipation, and do not decouple in the 
prograde direction. Requirement for dissipation 
underscores that a sufficient viscosity, and, there- 
fore, a presence of the gas is a necessary condition 
for prograde decoupling of nested bars. 

We note, that approximate constancy of fl s a 
for the decoupling bar discussed above has some 
analogy with the secular evolution of stellar bars. 
Those are known to extend to their corotation 
(e.g., Athanassoula 1992). Because the corota- 
tion radius typically lies in the flat part of the 
galactic rotation curve and so £1 ~ r _1 , this leads 
Q p a ~ const. One can show that the same relation 
can be worked out in the epicyclic approximation 
for decoupling gaseous bars as well, on the assump- 
tion that their size is limited by the radius of the 
IILR — which is supported by our simulations. 

An interesting question is to what degree both 
bars in the decoupled phase exchange energies 
via mechanism known as a mode coupling (e.g., 
Sygnet et al. 1988; Massett & Tagger 1997). The 
corotation radius corresponding to the final O s is 
positioned close enough to the OILR of the pri- 
mary bar, and their nonlinear interaction here is 
not out of question. However, mode coupling does 
not seem to contribute to the nested bar decou- 
pling itself. It can of course play an important 
role in locking these two modes in the resonance 
interaction when both bars are fully developed. 
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This issue will be addressed elsewhere. 

We note two more comments. First, the initial 
retrograde (in the frame of the primary bar) swing 
of the secondary bar does not seem to be of a prin- 
cipal importance. When the gas mass is increased 
just by 20%, the bar decouples directly into the 
progradc direction. Second, the model gaseous bar 
is expectedly triggered by the primary bar. In the 
control run with the same mass distribution but 
without the non-axisymmctric contribution from 
the large-scale bar, the nuclear bar did not develop 
and the central kpc remained nearly axisymmet- 
ric. 

In summary, in a fully developed galactic disk 
with a double ILR, a secondary gaseous bar forms 
in response to the gas inflow along the stellar 
bar. Initially corotating with the primary bar, this 
gaseous bar extends to radii where nuclear rings 
are known to form. Subsequent mass inflow and 
loss of circulation angular momentum to the stel- 
lar bar via gravitational torques strengthens this 
bar — a process which leads to development of 
chaotic orbits within the bar and, unlike in purely 
stellar bars, to a rapid gas infall of ~ 1 M Q yr -1 
over <~ 10 8 yrs across the central 200 pc, as gas 
cannot populate intersecting orbits. As the bar 
pattern speed appears to be inversely proportional 
to bar's shrinking size, the pattern speed increases 
manyfold, triggering a runaway dynamical decou- 
pling of the gaseous bar from the stellar bar. The 
gaseous bars extend only to their ILRs, unlike 
their stellar counterparts. Although numerical 
simulations of this avalanche inflow neccsarily are 
limited by a finite dynamic range which cuts off 
the inflow, it can in principle proceed to much 
smaller scales and has implications for fueling nu- 
clear starbursts and AGN. The outcome of this 
process depends on a number of physical processes 
which are not included in numerical simulations 
here, such as a runaway star formation and the 
degree of dumpiness of the interstellar matter. 
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